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The phenomenology of the scaling behavior of higher order structnre fnnctions of velocity differ¬ 
ences across a scale R in turbnlence should be built around the irreducible representations of the 
rotation symmetry group. Every irreducible representation is associated with a scalar function of 
R which may exhibit different scaling exponents. The common practice of using moments of longi¬ 
tudinal and transverse fluctuations mixes different scalar functions and therefore may mix different 
scaling exponents. It is shown explicitly how to extract pure scaling exponents for correlations 
functions of arbitrary orders. 
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Traditional measurements of anomalous scaling in tur¬ 
bulence are based on hot wire technology which yields 
information about the longitudinal components of the 
velocity field u{r,t) |p. Accordingly it is customary to 
consider the structure functions of longitudinal velocity 
differences: 

= , ( 1 ) 

Su{r, R,t) = u{r + R,t) — u{r,t) , (2) 

6 ue{r,R,t) = 6 u{r,R,t) ■ R/R . (3) 

It is well known that these structure functions appear 
to scale with scaling exponent which are anomalous 
(nonlinear functions of n): 

Sn{R) ~ R^- . (4) 


consider moments that are invariant under rotations [H; 
such invariants are expected to scale with pure scaling 
exponents that can be extracted from the data. 

The problem of mixing of different scalar functions 
does not exist for the second and third order moments of 
the longitudinal and transverse components. It is worth¬ 
while to go in detail through the analysis of the 2nd order 
moment in order to see why the longitudinal and trans¬ 
verse components are not a good choice, and why at the 
end it does not matter at this order. In an isotropic ho¬ 
mogeneous medium without helicity (with inversion sym¬ 
metry) the relevant symmetry group is the rotation group 
SO(3) whose irreducible representations can be expressed 
using the spherical harmonics Yi^rn- The most general 
form of the second order moment of velocity differences 
has contributions from £ = 0 and 2: 


Only recently has it become feasible, due to advances in 
experimental technology , and even more so in com¬ 
putational methods to measure other components 

of the velocity field. In particular a number of groups 
have focused on the transverse components 

6 ut{r, R, t) = 5u{r, R, t) — 5ui{r, R, t)R/R. (5) 

These groups studied the scaling exponents of the trans¬ 
verse structure functions 

T„(i?) = (|5wt(r,i^,^)r)-i^^" . (6) 

Two sets of measurements appear to imply that the scal¬ 
ing exponents ^ are the same as within experimental 
uncertainty J0,|| , whereas other numerical and ex¬ 
perimental studies indicate the opposite, i.e. that (f, 
are significantly smaller than for n > 4. 

The main point of this Letter is to demonstrate that 
higher order structure functions of longitudinal and 
transverse moments are not likely to exhibit clean scal¬ 
ing behavior, since they mix different scalar functions of 
R which may scale with different scaling exponents. In 
experimental and numerical studies in which all the com¬ 
ponents of the velocity field are available it is advisable to 


{5u°‘{r,R,t)5u^{r,R,t)) = 6 a/ 3 ao{R) 
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The coefficients in this expression carry the index £, mul¬ 
tiplying terms that are irreducible representations of the 
rotation group of dimension 2lA-\. The dimension of the 
irreducible representation is the number of tensor com¬ 
ponents that transform to one another upon rotation of 
the system of coordinates. All the tensor components of 
a given irreducible representation with a given value of i 
must have the same coefficient which depends only on R. 
On the other hand the scalar functions ao{R) and 02 ( 1 ?) 
may have different scaling exponents. 

Computing now the longitudinal and transverse mo¬ 
ments we find 

IDOL JD/S 

(SueSue) = {5u<^5u^) = oq (i?) - 203 (i?), ( 8 ) 

{5ut ■ 6ut) = (|(5mP) — {6ue6u() = 2ao{R) + 2a2{R). (9) 


Obviously these moments mix the two scalar functions 
with different weights. Fortunately the incompressibil¬ 
ity constraint forces ao{R) and 03 ( 1 ?) to have the same 
scaling exponent. We compute 
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meaning that the two functions must have the same R 
scaling, and therefore also the 2 nd order longitudinal and 
transverse components scale with the same exponents. 

The purity (and identity) of exponents of longitudinal 
and transverse fluctuations also holds for the third or¬ 
der moments. The most general form of the third order 
tensor {5u°‘5u^Su^) has contributions from ^ = \ and 3: 

{5u°^5u^5v?) = bi{R)[5o,pR? + 5^^R<^ + 

+b^{R)[5o.pR^ + 5c,-yR^ + 50^R°‘ - bR’^R^R^/R^] . 


{{5utf) = -^R'^R^R'^R^{5u‘^5u^5u'^5u^) (15) 

{\5ut\^) = {{5ut ■ 5utf) (16) 

DO: D/3 D7 D(5 

= - -^){5u‘^5u^5u^5u^) . 

A calculation yields 
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We again have two distinct scalar functions, each mul¬ 
tiplying a rotationally invariant form, and scaling with 
potentially different scaling exponents. Nevertheless, the 
incompressibility constraint provides one relation among 
the scalar functions, leaving us with one unknown. Kol¬ 
mogorov showed ||] that the rate of energy dissipation 
fixes the value of the remaining unknown. The form of 
the 3rd order tensor is thus fully determined, and a calcu¬ 
lation shows that (|(5Mf(r, i2, t)p) ~ (^\6ut{r, R,t)\^) ~ 
R. 

The first nontrivial example is the 4th order tensor 
{6u°‘6u^6u^Su^). The most general form of this tensor 
has contributions with £ = 0,2 and 4: 

{Su‘^Su^Su'>6u^) = co{R)Df^^ + c2{R)Df^^ 

+ Ci{R)Df^^ , (11) 


where 
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(13) 

(14) 


We see that in this case we have three independent scalar 
functions of R, i.e. co(i?), C 2 (i?) and ca{R), which in 
principle may have different scaling exponents. In this 
case the incompressibility constraint furnishes no rela¬ 
tion between these functions; the reason is that there ex¬ 
ist contributions in this tensor like {u°‘{r, t)u^{r, t)u'*{r+ 
R, t)u^{r -I- R, t)), and the divergence of such a contribu¬ 
tion (with summation on any tensor index) is not zero. 
In fact, incompressibility no longer places constraints for 
any of the higher order correlation functions for similar 
reasons. We note that there is no known way to justify 
why the three scalar functions should have the same de¬ 
pendence on R. We can compute now the longitudinal 
and transverse 4th order moments: 


We see that these components mix the three scalar func¬ 
tions with different coefficients. There are two possibil¬ 
ities: either all the scalar functions have the same lead¬ 
ing scaling exponent, or they have different scaling expo¬ 
nents. In the first case it is obvious that the longitudinal 
and transverse moments share the same scaling expo¬ 
nents. In the second case, for a sufficiently long inertial 
range, and for i? L where L is the outer scale of tur¬ 
bulence, the smallest exponent will dominate the scaling 
of both moments. Asymptotically the two moments are 
expected to have the same scaling behavior. However, 
if the three functions have different (leading) exponents, 
data with limited scaling range may lead to the erroneous 
conclusion that these moments have different scaling ex¬ 
ponents. It should be stressed that the amplitudes of the 
three scalar functions may be not universal, and different 
experiments may lead to different weights in this mixed 
representation. This may lead to a possible confusion or 
to conflicting results as seen in refs. [Ij- 0 . 

The more rational procedure that presents itself in 
light of this discussion is to compute the scaling behav¬ 
ior of the invariant scalar functions which are associ¬ 
ated with the higher order tensors. To achieve this we 
use the orthonormality of the irreducible representations, 
and observe that 


co{R) = Df^\5u°‘5u^5u^5u^) , (19) 

C2(i?) = Df^^du^^SuHuHu^) , (20) 

C4(i?) = Df^\5u°‘5u^5u'^5u^) . (21) 

Using the explicit form of the irreducible representations 
(|l^)-([l^) we can evaluate these functions and find 

co(i?) oc (|<5«|4Po ) OC (|5 m|4) , (22 ) 

a(S)oc(|H‘ft(^)) 

OC {\5u\^ [3{Suey - |i5mP]) , (23) 

C4(i?) oc {\6u\'^Pi 

oc {856ug — SOSujldul^ + 3|i5m|'^) , (24) 
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where Pt, are the standard Legendre polynomials of order 
We see that our scalar functions can be represented 
as particular combinations of transverse and longitudinal 
fluctuations. With data from a turbulent field it(r, t) one 
can compute in this way each of the independent scalar 
functions. Plotting them in double logarithmic plots (to 
get rid of the nonuniversal amplitudes) one has a good 
chance of extracting pure scaling behavior. After doing 
so one can return to the analysis of the longitudinal and 
transverse components with some understanding of the 
leading and subleading scaling exponents, to control the 
apparent scaling behavior in limited scaling ranges. 

These considerations are readily extended to higher or¬ 
der moments. The nth order tensor of velocity differences 
across a scalar R will have n/2 -|- 1 invariant scalar func¬ 
tions for n odd, and (n -|- l)/2 invariant functions for n 
even. There is no need to write down the explicit form 
of the irreducible representations, since the structure ex¬ 
hibited by Eqs.(p^ -(p4|) repeats at all orders. In other 
words, the independent scalar function d^{R) which is 
the function associated with the irreducible representa¬ 
tion of order ^ in the nth rank tensor of velocity differ¬ 
ences can be written in general as 

, ^<n, (25) 

where I has the same parity as n. Thus by simply exam¬ 
ining the Legendre polynomials in any textbook one can 
determine the precise combination of longitudinal and 
transverse fluctuations that is expected to scale with a 
pure exponents for any order n. 

In conclusion, it appears extremely worthwhile, in light 
of the growing abundance of high quality data on full tur¬ 
bulent velocity fields, to implement the approach detailed 
above. Since one confronts limited scaling ranges in most 


applications, it is mandatory to attempt to separate lead¬ 
ing from subleading scaling contributions in order to be 
able to make substantial conclusions about the numeri¬ 
cal values of scaling exponents. The procedure outlined 
above goes some way in this direction. 
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